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Abstract 

The optimal and minimal measuring strategy is obtained for a two-state system pre- 
pared in a mixed state with a probability given by any isotropic a priori distribution. We 
explicitly construct the specific optimal and minimal generalized measurements, which 
turn out to be independent of the a priori probability distribution, obtaining the best 
guesses for the unknown state as well as a closed expression for the maximal mean aver- 
aged fidelity. We do this for up to three copies of the unknown state in a way which leads 
to the generalization to any number of copies, which we then present and prove. 



*e-mail: guifre@ecm.ub.es 



1 Introduction 



A measurement allows us to extract only a small amount of the information needed to 
specify a quantum state. If our preparing device produces several identical copies of the 
unkown state, then measurements allow to extract more information, although only in the 
limit of infinitely many copies we acquire complete knowledge of the unknown quantum 
state. Performing an optimal measurement, the one which extracts the maximal possible 
amount of information about the state, and among these a minimal measurement, the 
one with the minimal number of outcomes, is always a priority, especially if the process 
leading to the state is rare or costly. It is also the broad subject of this paper. 

There are two aspects which significantly quantify the difficulty of the problem. One 
of them is the dimension of the Hilbert space which corresponds to the physical system we 
are considering. We will take the lowest one, two. The second is the a priori probability 
distribution function of the unkown state. If the state is known to be pure, the problem has 
been solved |]T], |^, |3| . The averaged, mean fidelity of the optimal measurements performed 
on N copies of a pure state is jl| 

FZipure) = £±1 (1.1) 

and the minimal measurements correspond, for N = 1 to 5, to || 

n^l(pure) =2,4,6,10,12 (1.2) 

outcomes. The aim of this paper is to solve this problem when we enlarge the a priori 
probability distribution function to include mixed states. More specifically, when one 
assumes that it is isotropic and otherwise arbitrary, but known. 

On the other hand, the difficult and heavily discussed issue about which is the abo- 
lutely unbiased probability distribution in the space of density matrices is not settled, 
and it might even not have an unbiased solution. In any case an unbiased distribution 
will be isotropic in the three-dimensional Poincare sphere covered by the Bloch vector 
which parametrizes the unknown density matrix and thus our results will be valid for any 
author's prefered candidate for an unbiased probability distribution. We will not discuss 
this issue further. 



Let us now outline the strategy defining optimal minimal measurements. We consider 
the simplest possible quantum system, a two-state system. It might be the spin of an 
electron, the polarization of a photon, an atom at very low temperatures so that only the 
two lowest hyperfine states matter, a linearly trapped ion for which only the ground and 
the first excited vibrational states are important, etc. This state is described by a 2 x 2 
density matrix 

p(b) = \ (I + b ■ a) = i±* | b)(b | | -b)(-b |, b =| b |< 1 , (1.3) 
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where b is the Bloch vector and | b) and | — b) are the eigenstates of p(b). These density 
matrices are prepared according to a known, isotropic, a priori probability distribution 
function given by 

/(6)>0 , 4tt j'dbb 2 f(b) = l. (1.4) 

Jo 

We will analyze the generalized measurements performed on the state corresponding 
to N copies of p(b), that is, p(b)® N , and determine which ones are optimal. There are 
two aspects to an optimal measurement: which are the positive operators correlated to 
the different outcomes, and which are the guesses which one makes, given an outcome, 
about the unkown state (which we shall call pi). Optimal measurements have to answer 
both questions by demanding that the guesses on average lead to the highest fidelity 
estimation of p(b), after averaging over the known probability distribution function f(b). 
We will then determine which of these optimal measurements are minimal, i.e. have 
the minimal number of outcomes. For more than one copy, N > 1, measurements may 
be collective and thus may involve entanglement. We will have something to say also 
about the relation between optimality and entanglement. The role of cloning as part 
of an optimal measurement will also be studied. We will also show that for more than 
two copies optimal measurements which are minimal are not complete, i.e. they involve 
positive operators with rank larger than one (and, yet, are optimal!). 

These are the main issues which will be presented for N — 1 to 3 copies in the next 
three sections. In sect. 5 we present and prove our general results for any N. The last 
section briefly recollects our findings and conclusions. 



2 N = 1 

Let us start with one single copy of p, N = 1, and use this example to present some of 
the systematics of our approach. 

We will first perform a generalized measurement [4] on p(b) with n outcomes, given 
by the operator sum decomposition 

X)4Ui = Ecfo = / , p i = p\>0 , Trpj = 1 (2.1) 
i=i i=i 

which implies 

n n 

Ec«=2 , E^1 = 0, (2.2) 

i=i i=i 

where is the Bloch vector of pi. If the outcome i is obtained, which happens with 
probability 

c*Tr(p(b) Pl ) =<*±(l + b-i£) , (2.3) 
one proposes pi as a guess for the unkown state p(b). The fidelity, i.e. the measure of the 
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goodness for a proposed guess, is quantified by 
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1 



F (p, p t ) = (Tr ^2^1/2 j =-(1 + 6.3 + Vl^V 1 " r * J . ( 2 - 4 ) 
where r« is the Bloch vector of pi. Thus, the fidelity averaged over all outcomes is 

F^\p) = i (l +b-s^j (l + b.f i + VrrpJY^) , (2.5) 

where the superscript reminds us that we are dealing with only one copy. From here 
the mean fidelity, i.e. the fidelity averaged over all unkown states p(b) weighed with the 
known probability distribution function f(b), is readily obtained 



T (N=1) = / dnf dbb 2 f(b)F (N=1 \ P ) 

= 7T dbb 2 f(b)Y,C 2 l [l + -S t 



VT^lPJl-r 2 ) . (2.6) 



i /l-6 2 \ a 
I a = 4vr / db b 2 f(b) — — , Jo = 1 , (2.7) 



With the notation 

F = 4vr f 1 dh b 2 f(b) (- 4 
(note that I a — 4/ a+1 > 0) the averaged fidelity reads 

F iN=1) = \ t c? (l + ^(1 - 4/i) * ■ n + 2J a/2 0^f) • (2.8) 

We have now to settle which is the best guess for the unknown initial state based 
on the result of our measurement, that is the proposed pi which leads to the highest 
mean fidelity. Let us first dispose of the case 4ii = 1, which corresponds only to /(&) = 
^p-lim e ^ ^(& — e), e > 0. It implies a vanishing Bloch vector and thus p{b) = \l, the 
completely random state. Since the unknown state is necessarily the completely random 
state, the state is known without performing any measurement whatsoever. We will thus 
always assume AI\ < 1, and only use 4Jx = 1 as a check-up of our results. Then from eq. 
( |2.8j ) maximization implies that the best guess corresponds to 

(1-4/05 (M) 



36/% + (1-4/0^? 



Notice that Pi + 1 Pi , but Pi is a known function of pi, as its coefficients depend only 
functionally on /(&). As f(b) is known, eq. ( |2.9|) determines the optimal guess in terms 
of pi. Substituting one obtains 



msvF > = FT' = 1^(1 + ^SQI 2 1/2 + (1 - 4/ 1 ) 2 ^ ) . (2.10) 
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We now have to determine the best measuring strategy, the one which leads to the largest 
possible fidelity. It is obviously given by s, = 1, i.e. by outcomes associated with rank-one 
projectors, and gives 

max^ =1) = F { Z: 1] = \ (l + \ ^QI 2 1/2 + (1 - 4A) 2 ) . (2.11) 

This is our result for one single copy of the physical system in state p(b) with a priori 
probability distribution fib). 

Notice that we have found that optimal measurements require necessarily an operator 
sum decomposition in terms of rank-one projectors. It is of course obvious that one can 
always perform an optimal measurement with rank-one projectors. Suppose, for instance, 
that we have some optimal operator sum decomposition with one operator of rank two, 
say pi. Then from its spectral decomposition 

Pi=Pi\pii)(pn\ + (l-Pi)\pi2)(pi2\ , (2.12) 

and from eq.( |2.3| ) 

c 2 Tr(p(&) Pj ) = c^Tr (p(b)\ Pil )(pii\) + c 2 (l - Pi )Tr (p(b)\ P v) (p a \) , (2.13) 

it is clear that taking as the guess for p for both outcomes associated to \pn) and \p i2 ) 
precisely pi, one can trade pi for its two rank-one eigenprojectors, having thus a measure- 
ment with only rank-one projectors. This result can be trivially generalized to iV copies 
and is of course well-kown [6]. We will use it without futher comments in obtaining F^ x , 
but it does not allow to analyze optimal measurements which are minimal, which will 
need a separate treatement. 

In the case we are considering here, N = 1, the outcomes are necessarily associated 
to rank-one operators and thus, from eq. ( |2.2|) , a minimal optimal measurement requires 
two outcomes, h.„j![~ = 2. This corresponds to a standard von Neumann measurement, 
which is a result unique for N — 1. For N > 1 optimal measurements are generalized 
measurements. 

A limit of interest corresponds to considering pure states, which is obtained by taking 
f(b) = j^liTn{, ^,x6(b — bo), b < 1. It follows that F^^ (pure) = |, which is the known 
result given in eq. fll.lj) . Notice that in this case pi = pi and thus the guess is precisely the 
pure state corresponding to the projector, while we have found that for mixed states the 
guess pi is a mixed state, different, though related, to the pure state corresponding to the 
projector. This is a new feature of optimal measurements. The two guesses correspond 
to two points in the interior of the Poincare sphere and symmetric with respect to its 
center. In the other extreme, discussed after eq. Q2.8[ ), when one knows that p(b) is the 
completely random state, we obtain F^^ 1 (random) = 1, as it should. One could think 

that minimizing F^^ with respect to f(b) would lead to |, as pure states cover the 
border of the Poincare sphere and thus maximize the naive distance between the states. 
This is not so. The probability distribution function f(b) = %h$(&(b) + 95(b — 1)) gives 

Fmax ^ = |(1 + ~7fo) < § an d we believe it to be the absolut minimum. 
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3 N = 2 

— * 

We will now study the situation in which two copies of the unknown state p(b) are avail- 
able, i.e. we have the state p(b) <8> pip). As we shall see, collective measurements appear 
here for the first time. 

Notice that defining the exchange operator V by 

V\<p)® \if>) = \<p), V = V* = V- 1 , (3.1) 

we have the following exchange invariance 

V(p® p)V = p® p . (3.2) 

We will consider generalized measurements for which outcomes correspond to rank-one 
projectors, as our purpose now is to build an optimal measurement. Thus the operator 
sum decomposition will be written as 



]Tc 2 |^>(^| = /, \iPi)eC 2 ®C 2 . (3.3) 

Given one decomposition one can obtain other decompositions as follows. First, obviously, 

n 

Ec? VhMtfclVW. (3.4) 



i=i 

Then, introducing the eigenstates of V built from \ipi) and V\ipi), 

\A)± = - I i== =(|^> ± V\A)) , (3.5) 

and, as 

iViX^i + vmmv = (i + ^m^)) \iPi) ++ m + a - mvm) m-m , (3.6) 

we have another decomposition 

Jt^^ + WW) \1>i)++<1*\ + (i - M |Vi>- -WW- (3-7) 
z i=l 

If the decomposition eq. ( |3.3| ) corresponds to an optimal measurement, so does eq. (|3.5|) 
just recalling eq. (|3.2|) and using the same guesses. Furthermore, as the probability of the 
i-th outcome is the sum of the probabilities of the i + and z_ outcomes of the decomposition 
ofeq. Q, 

c 2 

4 (ipi\p®p\if>i) = jO- + {<i>i\v\ii>i)) +<Vt|p 

c 2 

+ y(l- < ^l^lV'i >)-(Vi|/>®/#i>- (3.8) 
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it is enough to associate again the same guess to the i+ and i_ outcomes to make the mea- 
surement of eq. (|3.7|) optimal too. Thus optimal measurements can always be obtained 
by projecting on eigenstates of V. 

An equivalent way of presenting these results, and which will be more convenient for 
A^ > 2, is based on the identity 

V = S 2 - I (3.9) 
relating the exchange operator with the square of the total spin operator, 

S= -(a®I + I®a). (3.10) 

Eq. (|3.2|) now reads 

~S 2 ,p®p]=0 (3.11) 
and our previous results allow to write eq. ( |3.3| ) as 

n-l 

HH + E C ? \Ti)(Ti\=I (3.12) 

i=l 

where \a) is the singlet or antisymmetric state, and |rj) are triplet or symmetric states. 
This is an important result. It states that decomposing the Hilbert space of the two copies 
A and B into a direct sum of eigenspaces of S 2 , 

H (N=2) =u A ®H B = E @E 1 , (3. 13) 

where E s corresponds to the eigenvalue s(s + 1) of S 2 , it is enough to find optimal 
measurements in each of the spin eigenspaces for obtaining an optimal measurement in 
the whole space. The generalization of this result to A^ > 2 will be essential. It will then 
also be convenient to use both spin and exchange invariances simultaneously. 

We are ready to resume our general strategy for performing optimal measurements. 
First, the probability that the outcome corresponds to the singlet state is 

1-b 2 

(a\p®p\a) = — ^— . (3.14) 

For the triplet states we have found it convenient to use the Hilbert-Schmidt parametriza- 
tion 

|Tj) (ti| = - (l ® I + U ■ a ® I + / <g> ti ■ 3 + ii ■ a ® ii ■ a 

+ \Jl — t 2 (ui ■ a ® Ui ■ cr — %i ■ o ® %i ■ (?)^ (3.15) 

where ii, iii and Vi are n — 1 triads of orthonormalized vectors. Notice that ti is the Bloch 



vector of the reduced density matrix 

• / > I r M r = ^ 



Tt A \ri)(ri\ = Ti B \Ti)(Ti\ = - ii +ti- a) = pi , (3.16) 
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where we use subscripts A and B to earmark the Hilbert space over which the trace is 
performed. Furthermore from eq. (|3.12|) we have 



n— 1 n— 1 



£c 2 = 3 , £c 2 t1 = ( 3 - 17 ) 



i=i i=i 



and further restrictions on Vi and ij which will not be needed here. The probability 
that the outcome corresponds to \rj) is 

c? <T<|p ® Pki> = j (l + 2 6 • £ + (6 • U? + V 1 - ((& • «i) 2 - (6 • ^) 2 )) • (3-18) 

Once outcome i is obtained one proposes pi as a guess of the unknown state p(b) . From 
eq. ( |2.4j ) one obtains for the fidelity averaged over outcomes 

1 



F {N=2 \p) = -(l~b 2 )^l + b-f n + VT^¥^l-rl 

8 8=1 ^ 



l + b-fi + Vl-b 2 Jl-rf) . (3.19) 



The mean fidelity is obtained after averaging over the state space with the probability 
distribution function and reads 



F<"-> = i(/ 1+ 2/ 3/2v /I^i 

1 n_1 / 1 _ I \ 

+ gEc^l-/ 1 + -(l-4/ 1 )^-rl + 2(/ 1/2 -/3/ 2 ) V / l-^j • (3-20) 

From here the best guesses are readily obtained 

r n = (except for f(b) = -7-5(b — 1) when r n is not determined) 



16(/i/ 2 - / 3 / 2 ) 2 + (1 - 4/a) 2 t; 



=*i i = l,...,n-l (3.21) 



As before for N — 1, again p 4 ^ p^ is a function of p«, in fact a mixture of pi and the 
completely random state. Substituting the best guesses we obtain 

i i n ~i /i \ 

^T 2) = o J i + ^3/2 + «E c i + Tj^i/a-W + a-^i) 2 *?) • (3-22) 



6^ 1 V '2 



The best measurement strategy is obtained for U = 1, so that pi is a pure state and jr*) is 
a product state, without entanglement. This is a reasonable result, since p®p has neither 
entanglement nor classical correlations, so that it would be surprising that projecting on 
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entangled states would lead to an optimal measuring strategy. Notice also that this result 
of no entanglement, which we will reencounter later for N > 2, is independent of /(&). In 
fact, once the specification of the operator sum decomposition does not depend on /(&), 
it has to correspond to an optimal measurement strategy valid for pure states. But this 
is known [1, 2] to precisely require product states. For the singlet, which is a maximally 
entangled state, there are no alternatives, and thus the previous argument is irrelevant. 
The final result is 



F { l~ x 2) = \ + h/2 + yi6(h /2 -h/2) 2 + (l-^i) 2 (3.23) 

This final result reproduces the known limits. Indeed, the pure state result of eq. 
(|1.1| ) is readily obtained from eq. ( |3.23|) , when f(b) = -^b(p — 1). Also for the completely 



random state F^ ax (random) = 1. One can also check from the comparison of (F^ ax — |) 2 



for i — 1 and 2 that, as it should, 



-p(N=2) p(iV=l) . , 

max — max ■ yo.z,H:j 



Let us now analyze optimal measurements which are minimal. With the constraints we 
have been using for obtaining optimal measurements, i.e. an operator sum decomposition 
in terms of rank-one symmetric or antisymmetric projectors, the minimal n is five. This 
is because in the 3-dimensional symmetric (triplet) space a resolution of the identity in 
terms of symmetric product states needs four of them [3] , which together with the singlet 
makes five. When the unkown state is known to be pure, the outcome corresponding to 
the singlet never happens, and one can do with just four projectors. Let us now prove 
that one cannot do with less. 

Suppose we have an optimal measurement such that one of the rank-one projectors of 
its operator sum decomposition, with associated best guess p, is not symmetric nor 

antisymmetric. Obviously the best guess associated to VIV^^IV is also p. One can then 
build, following the arguments of eqs. (3.5-3~8|) an optimal measurement with 



and \ip)~ with associated best guesses p for both of them. But this is impossible, 
as we saw that the best guess associated to the antisymmetric state is the completely 
random state, while the one associated to the symmetric state has a non- vanishing Bloch 
vector (see eq. ( |3.21| ), and thus the best guesses cannot be equal. 

The very same reasoning forbids to have an optimal measurement with an opera- 
tor sum decomposition for which one of the operators has rank larger than one, as the 
associated rank-one projectors which appear in its spectral decomposition will have nec- 
essarily different best guesses. The upshot of all this is that for N = 2 minimal optimal 
measurements correspond to operator sum decompositions of rank-one symmetric or an- 
tisymmetric projectors, and thus have five outcomes, n^^f 1 = 5. We will see that for 
N > 2 the result that minimal measurements correspond to rank-one projectors does not 
hold. Notice that the five guesses are situated one at the center of the Poincare sphere 
and the other four on a concentric shell in its interior forming a regular tetrahedron. 

A related question to which we turn briefly is whether circumstances exist for which 
von Neumann measurements can be minimal and optimal. As C 2 <8> C 2 is of dimension 
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four a von Neumann measurement has four outcomes. We have seen that optimal mea- 
surements with four outcomes only exist when we know that the unknown state is pure. 
The question then is if the four triplet states, which are certainly not orthogonal, can be 
made orthogonal by adding them coherently to the singlet state. Notice that these states 
would not have a well-defined symmetry, but our previous proof that such states cannot 
be part of an optimal measurement fails precisely only for pure states, as then (cf. first 
of eq. (|3.21p ) r n is arbitrary. It is thus a legitimate question. Its answer is yes, for N = 2 
0. The answer for N > 2 is not known. 

Let us briefly go back to the situation in which we had one copy (section 2), and let 
us clone it with a state-independent universal quantum doner [7-11]. The conditions of 
strong [12] symmetry and isotropy of a universal l-to-2 quantum doner imply 

p(b) ^ pf ] = I + 7]{b-a®I + 1 ®b-a)+t ij a i ®a j ) 1 ty = t if - , (3.25) 

where r\ is the shrinking factor and where depends only on the vector b and the 
invariant tensor 5ij. Linearity, which originates in state- independence, and the absence 
of measurements in optimal cloning [13] forbids the quadratic dependence on 6j, so that 
eventually tij = t 5ij. It is also linearity which allows to clone straightforwardly for N — 1 
a mixed state by just mixing statistically the clones of the pure states which realize the 
mixed state. The values of the real parameters rj and t have to be such that p^ is a 
density matrix, i.e. such that its eigenvalues 

~(l±2br] + t) , \{l+t) , i(l-3f) (3.26) 

lie between and 1. Of course measuring on p^ will allow to learn the most about b for 
the largest rj possible. This is precisely what optimal cloning does: r] — | and thus t — g. 
We can now perform an optimal measurement on the optimal clone p^ , following closely 
the study of the N = 2 Vp^V = p^ . From the following results, 

Hp( 2 V> = , (r l \p^\r l ) = ^(l + b-t i ) , (3.27) 
the expression equivalent to eq. ( |3.19| ), after dropping an irrelevant part, is 



n— 1 



i=l 



F^(p) = -Y,cUl + b-q (l + 6.r1 + v / r^yi-rf) (3.2* 



This expression, together with eq. (|3.17|) , is identical to eq. (|2.5|), when eq. fl2.2| ) is 
recalled. We thus recover the result of eq. (|2.11|) . In words, optimal cloning can be part 
of an optimal measurement. As a byproduct we have checked that indeed p^ with t = | 
and T] = I is the optimal clone of p(b). 

Notice also the result shown in the first of eq. ( |3.27| ): the optimally cloned state 
lives in the triplet space. This is not surprising, as the singlet space cannot carry any 
information about the original cloned state. 
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4 N = 3 



Consider now three copies of the unknown state, p <S> p <8> p- Let us recall its exchange 
invariances 

[V A c,P®P®p] = Wbc,P®P®p] = 0, (4.1) 
where A,B,C are the subindices labeling the copies which are exchanged, and its spin 
invariances 

S 2 ,p®p®p] = [Sab,P®P®p] = °) ( 4 -2) 
where the partial and total spin operators are 

Sab = -(<t®I<8)I + I<8(t<8)I) , S = S AB + l -I ® I ® 3. (4.3) 
The first of eq. ( f4.2|) is obvious if one convinces oneself first that 

p®p®p = p 3 (S-b), (4.4) 



where Pn(x) is a polynomial in x of degree N. The second of eq. ([4.2|) follows then imme- 
diately. With the adequate generalizations in going from N = 2 to N = 3, it can be seen 
that in order to obtain optimal measurements it is enough to consider operator sum de- 
compositions whose elements are of rank one and project on states which are simultaneous 
eigenstates of S 2 and S\ B . Moreover these states should again be eigenstates of S ■ n for 
some n with maximal eigenvalue. Using the notation \s, sab, n), this leads immediately 
to the following states in terms of which the optimal operator sum decomposition can be 
built: 

3 

|-, l,n) = \n)\n)\n) 
|^,0, h) = \a)\n) 

\ll,fi) = -L(V A c-V BC )\<T)\fi). (4.5) 

The first state also corresponds to the completely symmetric representation of the permu- 
tation group generated by the exchange operators, and the other two correspond to the 
two-dimensional mixed symmetry representation of the same group. We may recall from 



ref. H that six states of the type of the first one of eq. (|4.5[) pointing into the six direc- 
tions of the vei 
maximal spin s 
decomposition 



tions of the vertices of a regular octahedron resolve the identity in the four-dimensional 

3 
2" 



maximal spin space, s — §. Therefore, we obtain the following optimal operator sum 



2 6 

(\ni)(ni\)® 3 + \a)(a\ <g> \n)(n\ + \a)(a\ ® | - n)(-n\ 



+\{V A c - V B c)\o){<t\ ® \n){n\{V AC - V BC ) 
+\{VAc-V B c)\<j)(<?\®\-h)(-h\{V A c-V BC ) = I. (4.6) 
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This result recalls the decomposition into eigenspaces E s<SAB oi S 2 and S AB , 

n (N=3) = H A ®H B ®H C = E ^ E ^ ^ 

and that under permutations Ei can be transformed into Ei v (Let us note here that the 
correctness of eq. ( |4.6|) has been confirmed by a brute force assumption-free computation 
which we performed in early stages of this work). Because of the isotropy of the probability 
distribution f(b) we just need to compute the following probabilities 



(h\(n\(h\p (g p (g p\n)\n)\n) = (h\p\h) 3 = - (l + b ■ 



n 



,2 



3 



i — b" 

{<y\{n\p (g) p ® p\o)\n) = (a\p (g p\a)(n\p\n) = — - — (l + b-nj 

~{cr\(n\{V AC - V BC )p ® p ® p{V AC - V B o)\a)\n) = {a\{n\p ® p ® p\a)\n) (4.8) 
where the last expression is obtained from 

\{Vac -V BC ) 2 \a)\h) = \a)\n) . (4.9) 
Putting all the pieces together we obtain (from eq. (|2.3[) ) 



F^ N=3 \p) = 1(1 - b 2 )(l + b ■ n) (l + b ■ f m + VT^Vy/l - rl 



+ -(1 + b ■ h) 3 fl + b ■ f s + VT 1 b 2 ^l - r 2 s j , (4.10) 

where r m and r s are the Bloch vectors of the proposed guesses of p corresponding to the 
mixed symmetry and completely symmetric projectors respectively. Angular integration 
over b leads to 

i^ =3) = l + l(/i -4/ 2 )n-F m + 2/3 /2 ^/l -r* 



+ (I 1/2 - 2/3/2) y/l-r* + —(3- Uh + 8/ 2 ) n-f s , (4.11) 



from which the optimal guesses are obtained for 

(/1-4/2) 



-n 



36 /2 /2 + (/l _ 4 j 2 )2 

3 - 14/i + 8/ 2 



/ 100(/ 1/2 - 2/ 3 /2) 2 + (3 - 14/i + 8/ 2 ) 2 
Substitution into eq. (|4. 10Q leads to our final result for N = 3, 

TdN=3) 1 1 



n . (4.12) 



+ — v /l00(/ 1/ 2-2/ 3/ 2) 2 + (3-14/ 1 + 8/2) 2 • (4.13) 
12 



This result reproduces the pure state result of eq. and gives 1 for the completely 

random state, as in previous cases. 

Let us finally come to those optimal measurements which are minimal. Up to now we 
have an optimal measurement with 10 outcomes. Remember that the only possibility of 
grouping together two rank-one projectors of the operator sum decomposition happens 
when the two different outcomes correspond to the same guess. Now from our results 
it is clear that this happens twice, that is the guesses corresponding to the 7-th and 9- 



th terms of eq. ( |4.6| ) are the same and given by the first of eq. ( 4.12 ), and the ones 
corresponding to the 8-th and 10-th terms of eq. Q4.6| ) are also the same and given by the 
first of eq. ( |4.12j ), but with opposite sign. Thus the minimal optimal measurement has 
eight outcomes, n^ in = 8. The corresponding positive operators On, s ,i and guesses Pn,s,i 
for iV = 3 are (cf. eq. (|4.fc>[) ) six for the space £"3/2,1: 



2 1 

03,3/2,; = ^\ni)(ni\ m , p 3 ,3/2,i = - {I + r s fn-a) , (4.14) 

and two for the space £1/2,0 © £1/2,1 : 

03,1/2,1 = \(T)(<T\®\n)(n\ + ^(V A c-V B c)\a)(a\®\h)(h\(V A c-V B c) 

P3,i/2,i = ^(I + r m n-a) 

03,1/2,2 = \o){<t\ <8> I - n)(-n\ + ^ (V AC - V BC ) W)(&\ ® | - h)(-h\ (V AC - V BC ) 

£3,1/2,2 = -(I-r m h-a) . (4.15) 

This is the first time in which a minimal optimal measurement has operators of rank 
two in its decomposition. The Bloch vectors of the corresponding guesses are situated on 
two concentric shells in the interior of the Poincare sphere. 

Notice that again the measuring strategy, i.e. eq. ( f4.6|) , is independent of f{b) and 
thus determined actually by what is known from [1, 2, 3]: for each s the pure state strategy 
for 2s copies is the optimal strategy. This will allow us to prove the general expression 
for an d n min f° r an Y N with relative ease in the next section. 



5 General results for N > 3 

We will analyze in this section optimal and minimal generalized measurements when a 
generic number N of copies of the unknown state are available. We present here the max- 

— (AO 

imal fidelity F max one can obtain on average by performing such collective measurements 
over p® N , together with the minimal number n^ n of outcomes an optimal generalized 
measurement can have. For any N we provide also a generalized measurement which is 
both optimal and minimal. Explicit results for the case N = 4 are worked out in order 
to illustrate the general expressions. 
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We first display our final, general results: 

= \ + E #^t^ f / ) Tj 9l (N,sy + g2 (N,sy, (5.i) 



where 



^i(jv, s ) = JdnJ Q dbb 2 f(b) j 



1 „ .o l-b 2 \ 2 1 + bA 



g 2 (N, s) = JdQ £ db b 2 f(b) (l-pj 2 " (i±^) " |, (5.2) 

6 2 is the third component of b and sq is (1/2) for even (odd) N. As for n^2 n we have 
found that 

N/2 



n 



E » (5-3) 



mm / j ' u ps 
s=s, 



where we define njjjp = n min (pure) , n^J = 1. For N = 1 to 5 this reads (using ||) 

^ = 2,5,8,15,20. (5.4) 

For N > 5 the minimal njjfp relies on a conjecture proposed in [3], and this is therefore 

also the case of n^ n for N > 5. 

For some very specific a priori probability distributions f(b) this number can be re- 
duced. This, though, corresponds only to cases in which there is an accidental degeneracy 
in the proposed guesses, as in the case f(b) = -^5{b — 1) (pure states). 

The optimal and minimal generalized measurements consists of the following decompo- 
sition of the identity operator in the space Jv- N > = C 2 ® N of the iV copies in terms of posi- 
tive operators On, s ,i and the corresponding guesses pn,s,%- for each se[s Q , s Q +l, y— 1, y], 
our optimal and minimal generalized measurement contains rS 2 ^ positive operators of the 
form 



if/Ll^i E V(\o)(o\^-°®\n s ,)(n s ,r 2s )^ , (5.5) 

1 \ 2 + S / iV • VeS N 



O - r 2 ^ 2S + 1 

2 t4tl \2 1 °/ ' VeS N 

where SV is the group of the N\ possible permutations of N elements acting on the Hilbert 
space of the N copies, and c 2 s i is such that 

N/2 nf a s) 

E E Nt . ti = I • (5.6) 

s=sq i=l 



The corresponding guesses are 



PN,s,i = ^(1 + r N , s n s>i ■ a), (5.7) 
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where 

92{N, s) 

rN,s = i (5.8) 
^ gi (N,sy + g 2 (N, S y 

The n£/' vectors Tig 2 dbTQ distributed according to their counterparts of the N = 2s 
case of optimal estimation of pure states as described in ||, and the coefficients c 2 s i satisfy 

(2s) (2s) 
Tips Hps 

£<i*W = , IX 4 = 2a + l. (5.9) 

i=l i=l 

For s = |, 1, |, | they are independent of i : c 2 ^ = ^fey- All these results are essentially 
unique. 

For N = A our results can be explicitly written as 



F { l x ' ] = \ + 2/ 5/2 + ^(2 - ll/i + 12/ 2 ) 2 + 36 (/ 1/2 - 3/ 3 /2 + h/2 

+ y(h-AI 2 f + 16 (13/2-/5/2) 2 (5.10) 

and 

^: 4) = 15. (5.11) 

The positive operator sum decomposition reads 

4 10 

/ = 04,0 + £ 04,1,* + £ ^4,2,* , (5.12) 
i=l i=l 

where to the rank- two projector 

o 4 ,o = 4 £ ^k)H ®k)Hvt (5.13) 
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there corresponds the guess 

1 



,0 = ^ Ko = 0). (5.14) 



The 4 rank-three positive operators 



04 lM = |j £ V|o-)<o-| ® |n M )<n M |® 2 V^" , 2 = 1,..., 4 (5.15) 

have associated guesses 

1 . /1 - 4/ 2 N 
P4,i,i = 77 (i + r 4 ,i • cr), r 4 1 = — == == (5.16) 

2 v /(/ 1 -4/ 2 )2 + 16(/ 3 / 2 -/5/2) 2 
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(here the are distributed according to a regular tetrahedron 0), and the 10 rank-one 
positive operators 



Oi,2,i = c^|n 2i )(n 2 . 



have associated guesses 

P4,2,i = + ?"4,2 "2,1 



^4,2 



i = 1,...,10 



(2 - 11/! + 12/ 2 ) 



(5.17) 



'(2 - ll/i + 12/ 2 ) 2 + 36(/ 1/2 - 3J 3 / 2 + / 5 /2) 2 

(5.18) 

(a concrete solution for n 2)i and c 2 j is given in ||). 

Let us now outline the proof of the above expressions. The proof will be based on a 
series of results which we have obtained along the previous sections and which we now 
put together in their generalized version: 

1. Permutation invariance 

For any element V of the permutation group of N elements, Sn, 



0, WeS 



N ■ 



(5.19) 



2. Spin invariance 

With the following notation for the composite Hilbert space, 

n (N) = n A ®n B ® ...n N , 

for the corresponding local spin operators, 

1 



Sa 
Sb 



A, 



-I® a® I 



®JV-2 



®)JV— 1 



s N = -I 

and for the partial and total spin operators 

M 

S m = J2S X , A<\JM<N 

x=A 



er, 



the spin invariances read 



S = S, 



(TV) 



0*2 n ® N ~\ — \ q 2 „® n 



2 „®jV 



^Ai P 



. 



(5.20) 



(5.21) 



(5.22) 



(5.23) 
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They are an immediate consequence of the following relatively straightforward result, 

p(bf N =p N (S-b) , (5.24) 

where Pn{x) is a polynominal of degree N in x. 

3. Direct sum decomposition 
Since 

r S 2 , Sf M) ] = \S? m ,Sf L) } =0 VM, L (5.25) 



the total Hilbert space can be written as a direct sum 

K {N) = ©.,{. (M) >E.,{. W > (5.26) 

where y are the eigenspaces of S 2 and Sf M ^,N > VM > A, with eigenvalues 

s(s + 1), {sm (sm + 1)} ordered with decreasing M, respectively. For instance, for N = 4, 

^i,f,i© E i,ii© E i,io© (s = l) 

£o,i,i ©^0,1,0 (s = 0). (5.27) 

Of course only those eigenvalues consistent with the spin composition rules appear. 
4. Permutation group equivalence 

For a given s < y all the spaces # s ,{s (M )} corresponding to it can be obtained from 
one of them with the help of the elements of the permutation group. The one which we 
retain for our proof as reference space is the one with the maximal number of vanishing 
partial spins, 

N 

# S)S -i, s -i,...o,i,o ( with y - s zeros )- ( 5 - 28 ) 

There are as many of these equivalent spaces as the dimension of the irreducible repre- 
sentation of Sn in a space of total spin s. 

'»<•> = ( f + J S^TTT < 5 - 29 > 

One can check the dimensional consistency of the previous expression from eq. ( 5.26|) , 

JV 

2 N = (2s + 1) d N (s) s = or - (5.30) 

S = SQ 



5. Optimal pure state measuring strategy 
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In each of the reference spaces of the type of eq. ( |5.28|) where any vector is of the form 

|cr>®£ - s <g> |V) , |V) e C 202s (5.31) 

the best measuring strategy turns out to be the one corresponding to 2s copies of an 
unknown pure state [1, 2, 3], and thus projects onto states of the form 

\af%- s <g> |n)® 2s (5.32) 

Notice that the singlets act as an identity in the reference space of eq. ( |5.28] ) and that the 
states flOp are the ones in eq. (|5.31| ) with less entanglement. From here, and recalling 



eq. ( |5.29| ), one readily obtains eqs. (|5.5| ) and (|5.6|). The fact that the guesses of eq. ( |5TT| ) 
can be grouped together due to the permutation equivalence, and thus have to be made 
only for the reference space, has been taken into account already in writing eq. ( |5.5D . 
Notice that the operators of eq. (|5.5|) are of rank cIn(s). 



We are now ready to perform the final computation of 

JV (2s) 

F^L =EE/*( db b 2 f(b)Tr (0 N , s>i p® N ) F(p, p N ^). (5.33) 



From 



l-6 2 \ 2 (1 + b-n. 



n(0 N ^°»)=<* t Ms)( L 1 L ) ^p^] , (5.34) 



which is obtained from eq. ( O ), Q5.5| ) and ( |5.29| ), eq. ( |5.33| ) can be written as 



A ,„ ^ , , s f ,~ f 1 „ ,9 ,„x /I - & \ 2 / 1 + 6-n 



^ = E( 2s + / d« / db b 2 f(b) 

S — So 

l - (l + 6 ■ h + v/I^fe 2 y/l - r\ 8 ) (5.35) 



where we have used eq. (|5.9|) , as the contributions corresponding to different i are the 
same, eq. ( |2.4| ) for the fidelity, and where the subindices of h Sii have been dropped, given 
their irrelevance at this stage of the computation. In eq. (|5.35 ) the first term gives ~ 



and the other two depend on rjv jS , which is fixed by maximization. Choosing h in the 
direction of the z-axis, and with the definitions of eq. (|5.2|) , one immediately obtains eq. 
( |5.S| ) and finally our main result eq. ( |5.1| ). The result refering to the number of outputs 



of minimal measurements, eq. (|5.3| ), follows from our point 5 above. 



6 Conclusions 

We have built the optimal and minimal measuring strategy for N copies of an unknown 
mixed state prepared according to a known, isotropic, but otherwise arbitrary probability 
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distribution. The strategy is universal, i.e. independent of the probability distribution. 
Except for one single copy, optimal measurements have to be generalized measurements. 
We have obtained a closed expression for the maximal averaged mean fidelity, and the 
associated minimal number of outcomes. In obtaining these expressions some interesting 
windfall results emerged. They are: 

1) Best guesses are not universal. They are pure states only if the unknown state is 
known to be pure. 

2) Optimal measurements require projecting onto total spin eigenspaces, and within 
each such subspace, onto total spin eigenstates with maximal total spin component in 
some direction. This allows to relate them with optimal measurements corresponding to 
a smaller number of copies of unknown pure states. 

3) Optimal measurements which are minimal have, beyond two copies, outcomes as- 
sociated with positive operators of rank larger than one and, beyond three copies, less 
outcomes than dimensions of the Hilbert space. These optimal measurements are thus 
incomplete! Completing them is useless. 

Our results also set the limits to optimal cloning of mixed states. The techniques 
developed here for dealing with copies of mixed states will be useful for solving related 
problems. 

After finishing this work we learned from Ignacio Cirac that he has done, together 
with Artur Ekert and Chiara Macchiavello, somewhat similar work using basically the 
same techniques. 
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